Transition from a strongly interacting ID superfluid to a Mott insulator 
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We study one-dimensional trapped Bose gases in the strongly interacting regime. The systems 
are created in an optical lattice and are subject to a longitudinal periodic potential. Bragg spec- 
troscopy enables us to investigate the excitation spectrum of the one- dimensional gas in different 
regimes. In the superfluid phase a broad continuum of excitations is observed which calls for an 
interpretation beyond the Bogoliubov spectrum taking into account the effect of quantum deple- 
tion. In the Mott insulating phase a discrete spectrum is measured. The excitation spectra of both 
phases are compared to the three-dimensional situation and to the crossover regime from one to 
three dimensions. The coherence length and coherent fraction of the gas in all configurations are 
measured quantitatively. We observe signatures for increased fluctuations which are characteristic 
for ID systems. Furthermore, ceasing collective oscillations near the transition to the Mott insulator 
phase are found. 

PACS numbers: 05.30.Jp, 03.75.Kk, 03.75.Lm, 73.43. Nq 



Quantum gases trapped in the periodic potential of an 
optical lattice have opened a new experimental window 
on many-particle quantum physics. The recent observa- 
tion of the quantum phase transition from a superfluid 
to a Mott insulating phase in a Bose gas T] has offered 
a first glimpse on the physics which is now becoming ex- 
perimentally accessible. However, the full wealth of pos- 
sibilities has yet to be explored. Besides controlling the 
effect of interactions in the trapped gas, it is conceivable 
to induce disorder, to change the dimensionality of the 
system, or to trap Fermi gases or Bose-Fermi mixtures. 
The realization of these systems is expected to provide a 
deeper understanding of general concepts related to su- 
perfluidity and superconductivity. 

Here we use the optical lattice to realize a strongly in- 
teracting Bose gas in one spatial dimension and to study 
the crossover to three dimensions. Emphasis is put on 
the measurement of excitation spectra which character- 
ize the transition from the superfluid 0, Q to the Mott 
insulating state P, 0, 01 • Several features observed in the 
spectra go beyond the description of current theoretical 
models. 

Degenerate Bose gases trapped in the lowest band of an 
optical lattice can be modelled using the Bose-Hubbard 
Hamiltonian 0, 0, 1 which the hopping of atoms 
between neighboring lattice sites is characterized by the 
tunnelling matrix element J, while the interaction en- 
ergy for two atoms occupying the same site is given by 
U . The physics of this model is governed by the ratio 
between U and J, i.e. between interaction and kinetic 
energy. This parameter can be controlled by changing 
the depth of the lattice potential. If the ratio U / J is be- 
low a critical value the atoms are superfluid. Above this 
value the system becomes Mott insulating. We access the 
one-dimensional regime 0, Hol ll]| using an anisotropic 
optical lattice consisting of three mutually perpendicu- 
lar standing waves. By choosing large potential depths 
in two axes we can selectively suppress tunnelling and 



hopping is possible only along one dimension, so that an 
array of one-dimensional tubes with periodic modulation 
along their axis is formed. 

We produce almost pure Bose-Einstein condensates of 
typically 1.5 x 10^ ^"^Rb atoms in the |F = 2,tof = 2) 
hyperfine ground state in a magnetic trap with trapping 
frequencies of = 27rxl8Hz, tUy = 27rx20Hz, and a;^ = 
27r X 22 Hz. The optical lattice is formed by three retro- 
reflected laser beams which are derived from laser diodes 
at a wavelength of A = 826 nm the position of 

the condensate the beams are circularly focused to 1/ e^- 
radii of 120 /um (x and y axes) and 105 /im (z). The 
three beams possess mutually orthogonal polarizations 
and their frequencies are offset with respect to each other 
by several ten MHz. The linewidth of the lasers is of the 
order of 10 kHz. 

In order to load the condensate into the ground state of 
the optical lattice, the intensities of the lasers are slowly 
increased to their final values using an exponential ramp 
with a time constant of 25 ms and a duration of 100 ms. 
The resulting optical potential depths Vx^y^z are propor- 
tional to the laser intensities and are conveniently ex- 
pressed in terms of the recoil energy — with 
k — ^ and the atomic mass m. To prepare an array 
of one-dimensional tubes, two lattice axes are ramped to 
V± = Vx — Vz ^ SQEfj and the third one to a much 
lower value Vaxfi = Vy. In this configuration the trans- 
verse tunnelling rates Jx and Jz are small and contribute 
a correction of the order of Jx,z/fJ- ^ 1 to the ID charac- 
teristics of the individual tubes, where fi is the chemical 
potential of the sample. It is convenient to include the 
anisotropic tunnelling between all nearest neighbor sites 
which yields the resulting J = 2{Jx + Jy + Jz)- 

We study the excitation spectrum by employing am- 
plitude modulation of the axial lattice potential Vax 
to perform two-photon Bragg spectroscopy [l^. The 
lattice potential takes the form Vaxiy,t) — {Vax,o + 
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-^mod Sin (27ri/mo(jt)) sin^(fcy). The modulation with am- 
phtude Amod and frequency i^mod introduces two side- 
bands with frequencies ±Vmod relative to the lattice laser 
frequency which define the energy hvmod of the excita- 
tion. Due to the Bragg condition atoms scattering two 
photons receive a momentum transfer of Ohk or 2hk. In 
contrast to applying a potential gradient across the lat- 
tice Q , this method is not susceptible to effects like Bloch 
oscillations and Zener tunnelling which occur for low ax- 
ial lattice depths. Furthermore, the excitation energy is 
precisely determined and does not involve any parame- 
ters that need calibration. 

After the excitation, the experimental sequence con- 
tinues by ramping down the lattice potentials linearly in 
15 ms to Vax — V± — 4: En where the atoms are able to 
tunnel again in all three dimensions between the sites of 
the lattice. To allow for re-thermalization of the system, 
the atoms are kept at this lattice depth for 5 ms. Then 
all optical and magnetic potentials are suddenly switched 
off. The resulting matter wave interference pattern is 
detected by absorption imaging after 25 ms of ballistic 
expansion. The width of the central momentum peak 
is taken as a measure of how much energy has been de- 
posited in the sample by the excitation. If the energy 
increase is small, the peak is well fitted by a bimodal dis- 
tribution. For resonant excitation there is only a single 
gaussian component, reflecting that the temperature of 
the atoms has significantly increased. To be independent 
of the shape of the peak we use the full width at half max- 
imum (FWHM) as a measure of the introduced energy. 
Although this underestimates small energy increases, the 
important resonances and features of the spectra are well 
shown. 

The duration tmod = 30 ms and amplitude A^od = 
0.2Vax,o of the modulation are chosen such that the re- 
sulting excitation of the condensate does not exhibit sat- 
uration effects for all measurements presented here. We 
have verified that all atoms remain in the lowest Bloch 
band by adiabatically switching off the lattice potentials 
[l^ after the modulation. When we load a cold thermal 
cloud into the lowest Bloch band and apply our modula- 
tion scheme, we do not observe excitations. 

Figure^displays the fundamental change in the excita- 
tion spectrum for a ID Bose gas when the crossover from 
the superfluid to the Mott insulating phase occurs: The 
broad continuum of the superfluid contrasts with the dis- 
crete spectrum of the Mott insulator. One surprising fea- 
ture is that we can excite the superfluid with our scheme 
at large hvmod contrary to predictions for the weakly 
interacting superfluid in an optical lattice formed by a 
single standing wave . In our experiment strong inter- 
actions lead to a significant quantum depletion, which is 
« 50% for the ID configuration even with U j J — 2.3 J^J. 
Therefore this parameter may not be regarded small as 
in standard Bogoliubov theory, but hig her order excita- 
tions should be taken into account 'iq. 
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FIG. 1: Spectroscopy of the ID superfluid (open circles) and 
the Mott insulating phase (filled circles) with values of [// J 
of approximately 2.3 and 14 respectively. Error bars reflect 
the statistical error of 5 measurements. 



In combination 



with the broken translational invariance in the inhomo- 
geneous trap, this could explain the non-vanishing ex- 
citation pro bability observed in the experiment at high 
energies |l7j . 

A full series of spectra for different values of U/J rang- 
ing from the superfluid via the crossover region to the 
Mott insulating phase is shown in figure[21 FigurelSt dis- 
plays the ID situation with V± = 30 Er. The crossover 
from one to three dimensions (figure ^p) is achieved by 
reducing the transverse confinement of the optical lattice 
to V± — 20E]i. Then the tunnelling time between the 
ID tubes is of the order of tmod and thus the interac- 
tion between the tubes is not negligible any more p^ . 
Finally, figure [3; shows the three-dimensional case with 

In one dimension (figurel^L) we observe the appearance 
of the discrete structure, which is characteristic for the 
Mott insulating phase, between U/J = 4 and U/J = 8. 
Above U/J = 20 there is no more background due to the 
vanishing superfluid component. Our results are in ac- 
cordance with the prediction U/J = 5.8 for entering the 
n = 1 Mott insulator based on a mean- field theory jQ, . 
Calculations beyond the mean-field approach give an on- 
set of the Mott insulating phase in the homogeneous ID 
system at U/J~ 1.8 However, the finite size of the 
trap prohibits a sharp transition [lOj , so that the fraction 
of Mott insulating atoms increases gradually with U/ J. 

For the superfluid we obtain spectra which differ sig- 
nificantly from the results of ref. [1| , since the superfluid 
excitations decrease at higher energies. This decrease is 
rather slow for the ID gas but becomes more pronounced 
when the tunnelling between the ID gases is increased 
(Figure l^h and Figure 13;). Our excitation scheme does 
not induce dephasing that occurs when the strongly in- 
teracting condensate is accelerated near the edge of the 
Brillouin zone |l9l | . This might cause the broadening and 
the background in the tilted lattice experiments at high 
energies in ref. 0. The width of the superfluid spectra 
for the ID gas is on the same order as twice the width of 
the lowest band for Bogoliubov excitations [20j . 
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FIG. 2: The measured excitation spectrum of an array of ID 
gases {V± — 30 Er) is shown in (a) for different values of 
Vax,o- The interaction ratios U/J given in brackets are calcu- 
lated numerically using a band structure model in the tight- 
binding approximation Q . Spectrum (c) shows the superfluid 
to Mott insulator transition in the 3D case {V± = Vax,o)- The 
crossover region between the one- and the three-dimensional 
system {V± — 20 Er) is shown in (b). 
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In the Mott insulating phase we find the first resonant 
peak for all data sets close to the calculated value of U. 
A second peak appears at (1.91 ± 0.04) times the energy 
of the first resonance, somewhat smaller than the value of 
2 reported in This resonance might be attributed to 
defects where lattice sites with n = 1 atom next to sites 
with n — 2 atoms are being excited. For the ID system 
and in the dimensional crossover regime (Fig. |3i and b) a 
much weaker resonance appears at (2.60±0.05) times the 
energy of the first resonance which could indicate higher 
order processes of two atoms tunnelling simultaneously. 
In figure Et- we plot the rms width of the first resonance 
in the excitation spectrum of the Mott insulating phase 
when fitted by a gaussian. In figure ^jp we show the 
ratio of the amplitudes of the second and the first peak. 
Apparently, in the ID system the first peak is wider and 
the second peak more pronounced as compared to the 
3D situation, which could be an indication of increased 
fluctuations in ID systems. 
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FIG. 3: (a) Width of the first resonance peak in the spectrum 
of the Mott insulator, (b) Ratio of the amplitudes of the sec- 
ond and the first peak of the spectrum of the Mott insulator. 
The error bars mark the error of the gaussian fits. 

Compared to the superfluid properties the coherence 
properties of the system provide complementary infor- 
mation about the state of the gas. They are probed by 
studying the matter wave interference pattern P, l2ll| . 
Here we first prepare the array of ID systems as above 
but do not apply our excitation scheme. Instead, af- 
ter holding the atoms at the final lattice depth for th = 
30 ms, we increase Vax rapidly (< 40 ^s) to about 25 Eji 
and then abruptly switch off all optical and magnetic 
trapping potentials. This procedure projects the differ- 
ent initial configurations onto the same Bloch state. To 
extract the number of coherent atoms Ncoh from the in- 
terference pattern, the peaks j23| at Ohk, ±2hk and ±ATik 
are fitted by gaussians (Fig. ^Jj). Incoherent atoms give 
rise to a broad gaussian background which dominates for 
higher Vax ft- Taking this fit as a measure of the number 
of incoherent atoms Nincoh, we calculate the coherent 
fraction fc = jr — . As shown in figure^, fc de- 



creases slowly to zero for increasing values of U/J and 
appears almost independent of the dimensionality. This 
coincides with the prediction that for strongly interact- 
ing Bose gases in optical lattices the superfiuid fraction 
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can be significantly different from the coherent fraction, 
and that the decrease of fc is not a sufficient signature 
of entering the Mott insulating phase In figure 0]: 

we plot the width of the central peak of the interference 
pattern, which is a measure of the coherence length of 
the gas. An increasing width is a good indicator for the 
presence of a Mott insulating phase since even a small 
Mott insulating domain reduces the coherence length of 
the sample, as elucidated in numerical calculations [23 |. 
Our data shows that the increase in width starts at much 
lower values of U/ J for the ID gas than for the 3D gas. 
This supports the expectation that due to the more pro- 
nounced quantum fluctuations in the ID geometry the 
gas enters the Mott insulating state at lower values of 
U/J HElIillil- Experimentally, thermal fluctuations 
may also contribute to the observed width. However, by 
comparing the width of the interference peak for differ- 
ent hold times {th = 1ms and = 30 ms) wc find that 
the primary effect of the additional heating is an overall 
increase of the width rather than a change of the slope 
of the curve. 
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FIG. 4: (a) Coherent fraction vs. U/J. The error bars are 
determined by the statistical error of 4 measurements, (b) 
The column sum of the optical density (circles) and the fits 
(solid line) from which the number of coherent and incoher- 
ent atoms are deduced for the ID case. The inset shows 
the absorption image after 10 ms of time-of-flight (dimensions 
360 fim X 467 urn) . (c) Width of the central momentum peak 
vs. U/J. 



One feature of the Mott insulating phase is its lack 
of compressibility 0. If a potential gradient is appHed 
across the lattice, the inability to redistribute atoms in 
the lattice results in a breakdown of the flow of atoms. 
Even if only a fraction of the atoms in our inhomogeneous 
trap is locally incompressible and pinned in the Mott in- 
sulating state mi, collective oscillations should suspend. 
In the experiment, the array of one-dimensional gases is 
prepared in the same way as before. Then, we apply 
a small potential gradient (far less than any resonance 



gradients 0) for 1ms to accelerate the atoms in the ax- 
ial direction. We observe a very strong increase in the 
damping of the small amplitude dipole oscillations even 
for very low Vax,o, a behaviour which is different from 
the three-dimensional gas. Critical damping is reached 
at Vaxfl — 3i?fl, and at Vax,o = lOSij we can not dis- 
place the atom cloud any more. 

In conclusion, we have prepared an array of one- 
dimensional strongly interacting Bose gases and mea- 
sured the excitation spectra in the superfluid and the 
Mott insulating regime. The spectrum of the superfluid 
exhibits excitations at high energies which are not pre- 
dicted by current theories for weakly interacting bosons 
in an optical lattice and have so far only been discussed 
in the context of superfluid Helium. This shows that the 
strong interactions and the significant quantum deple- 
tion change the properties of the Bose gas considerably, 
already much below the crossover to the Mott insula- 
tor. The effect of reduced dimensionality on the crossover 
from superfluid to Mott insulator was most pronounced 
in the measurements of the coherence length of the gas. 

We would like to thank H. P. Biichler and G. Blatter 
for insightful discussions, and SNF and QSIT for funding. 
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